Contact Equivalence Problem for Nonlinear Wave 
Equations 



Oleg I. Morozov 

Department of Mathematics, Snezhinsk Physical and Technical Academy, 

Snczhinsk, 456776, Russia 

morozov@sfti.snz.ru 

Abstract. The moving coframe method is applied to solve the local equivalence 
problem for the class of nonlinear wave equations in two independent variables under 
an action of the pseudo-group of contact transformations. The structure equations 
and the complete sets of differential invariants for symmetry groups are found. The 
solution of the equivalence problem is given in terms of these invariants. 

AMS classification scheme numbers: 58H05, 58J70, 35A30 
Introduction 

In this article we consider a local equivalence problem for the class of nonlinear second 
order wave equations 

w tt = f{x,w x ) w xx + g(x,w x ) (1) 

under a contact transformation pseudo-group. Two equations are said to be equiva- 
lent if there exists a contact transformation mapping one equation to the other. Elie 
Cartan developed a general method for solving equivalence problems for submanifolds 
under an action of a Lie pseudo-group, jT] - [3]. The method provides an effective 
means of computing complete systems of differential invariants and associated invariant 
differential operators. The necessary and sufficient condition for equivalence of two 
submanifolds is formulated in terms of the differential invariants. The invariants 
parameterize the classifying manifold associated with given submanifolds. Cartan's 
solution to the equivalence problem states that two submanifolds are (locally) equivalent 
if and only if their classifying manifolds (locally) overlap. The symmetry classification 
problem for classes of differential equations is closely related to the problem of local 
equivalence: symmetry groups and their Lie algebras of two equations are necessarily 
isomorphic if these equations are equivalent, while the converse statement is not true in 
general. The preliminary symmetry group classification for the class (JIJ is given in j^j. 
In [10] , it was proposed to transform equation (JTJ) to the equivalent quasi-linear system 
of the first order 



u t = a(x,u)v x , v t = b(x,u)u x , 



(2) 
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and the symmetry classification for non-linearizable cases of this system is given. In 
[To] several cases of infinite symmetry algebras for equation (0) are found, and one 
linearizable case is given. 

In the present paper, we apply Cartan's equivalence method, jTj - j^j, [E], [IS], in 
its form developed by Fels and Olver, [SI E|, to find all differential invariants of sym- 
metry groups and to solve the local contact equivalence problem for equations from 
the class (0) in terms of their coefficients. Unlike Lie's infinitesimal method, Cartan's 
approach allows us to find differential invariants and invariant differential operators 
without analyzing over-determined systems of PDEs at all, and requires differentiation 
and linear algebra operations only. 

The paper is organized as follows. In Section 1, we begin with some notation, and 
briefly describe the approach to computing symmetry groups of differential equations 
via the moving coframe method of In Section 2, the method is applied to the class 
of nonlinear wave equations ©• Finally, we make some concluding remarks. 

1. Pseudo-group of contact transformations and symmetries of differential 
equations 

In this paper, all considerations are of local nature, and all mappings are real analytic. 
Suppose £ = R n x M. m — > R n is a trivial bundle with the local base coordinates 
(x 1 , x n ) and the local fibre coordinates {u l , ...,u m ); then by J 1 ^) denote the 
bundle of the first-order jets of sections of S, with the local coordinates 
i G {1, ■■■,n}, a G {1, ...,m}. For every local section (x l , f a {x)) of S, the corresponding 
1-jet (x\ f a (x), df a (x)/dx l ) is denoted by ji(f). A differential 1-form i? on J x {£) is 
called a contact form, if it is annihilated by all 1-jets of local sections: = 0. 

In the local coordinates every contact 1-form is a linear combination of the forms 
d a = du a — pfdx 1 , a G {!,..., m} (here and later we use the Einstein summation 
convention, so pf dx % = Yh=i P? \ etc.) A local diffeomorphism 

A : J\E) -> J 1 ^), A : {x\u a ^) » (x*,^,^), (3) 

is called a contact transformation, if for every contact 1-form the form A*$ is also 
contact, in other words, if A*$ a = du a —pfdx 1 = (g(x,u,p) i? 10 for some functions 
on J^S). 

Cartan's method of equivalence, El EH], allows us to compute invariant 1-forms 
which define the pseudo-group of contact transformations. The result of its application 
is the following (see Consider the lifted coframe 

6° = a° (du p -p^dx j ), 

H* = 4 e p + b)dx j , (4) 
S? = ftp + Sg ^ + 4 Bl dp] 
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on J 1 (£) x TC, where Ti is the Lie group of block lower triangular matrices of the form 
/ a$ \ 



c*<# b) 



V (/£ + 9?k 4) 4 9% h) a% B{ I 
and the parameters a^, bj, c^, f^, and gfj obey the requirements det {a^j ^ 0, 
det (b l A 7^ 0, b\B k - = 5], and = g" { . Then a transformation T : ^(S) x Ti — > 
J 1 (£^) x 7i satisfies the conditions 

T* 0° = e a , T* S* = 5*, T* = £f 

if and only if it is projectable on J 1 ^) and its projection A : J l {£) — > J l {£) is a contact 
transformation. The lifted coframe has the structure equations 

do a = $° a e' 3 + E k A E£, 

ds* = #i a s fe + n 7 a e 7 , (5) 
dE? = $ 7 a ej — aEj + a$ a + n« a 

where 11^, A°g, and £H are 1-forms on J 1 ^) x and, as it is shown in [TT] . 

the coframe is involutive. 

The structure equations (0) remain unchanged if we make the following change of 
the modified Maurer - Cartan forms ^ k , IT , A^, and flfj : 

n; ^ n 7 + M* 7 ~ fc + iv 7e e e , (6) 

where K" L z kj , ML,, , P^ 7 , and are arbitrary functions on J 1 (£^) x Ti 

satisfying the following symmetry conditions: 

K" e = K^, L kj = L jk , N ye = N l n , 

pa pa f)Ot r\a pa pa pa 

r ifrl ~ r i,fii ^if3k ~ Wkfiii ^ijk ~ ^ikj ~ ^jik- 

Another approach to construct 1-forms characterizing contact transformations is 
presented in jHj. 

Suppose TZ is a first-order differential equation in m dependent and n independent 
variables. We consider TZ as a sub-bundle in J l {£). Let ContiTZ) be the group of 
contact symmetries for TZ. It consists of all the contact transformations on J 1 ^) 
mapping TZ to itself. The moving coframe method, OIZI, is applicable to find invariant 
1-forms characterizing Cont(7Z) is the same way, as the lifted coframe (gj) to J l {S) x Ti 
characterizes Cont{J 1 {8)). We briefly outline this approach. 

Let i : TZ —>■ J x {£) be an embedding. The invariant 1-forms of Cont(7Z) are re- 
strictions of the coframe (p£J) on TZ: 9 a = i*O a , £ l = t*E\ and erf = (for brevity 
we identify the map t x id : TZ x Ti -> J 1 ^) x Ti with l : TZ ^ J 1 (£)) . The forms 
£*, and erf have some linear dependencies, i.e., there exists a non-trivial set of functions 
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S a , Tj, and U % a on 1Z x TC such that 6* a + Tj + U % a erf" = 0. These functions are lifted 
invariants of Cont(lZ). Setting them equal to appropriate constants allows us to specify 
some parameters a^, 6j, cjg, and of the group TC as functions of the coordinates 
on 1Z and the other group parameters. 

After these normalizations, some restrictions of the forms 4>p = L*Qp, ip\ = 
7Tp = t*ITg, Xfp = £*A^g, and ufj = i*Qfj, or some their linear combinations, become 
semi-basic, i.e., they do not include the differentials of the parameters of TC. From ©, 
we have the following statements: (i) if 0^ is semi-basic, then its coefficients at crj and 
£ J are lifted invariants of Cont(lZ); (ii) if ip l k or it^ are semi-basic, then their coefficients 
at aj are lifted invariants of Cont(lZ). Setting these invariants equal to some constants, 
we get specifications of some more parameters of TC as functions of the coordinates on 
TZ and the other group parameters. 

More lifted invariants can appear as essential torsion coefficients in the reduced 
structure equations 

d9 a = <$>% A^ + £ fc A<r£, 

da* = <P«A a] - $ A at + A$ A ^ + w?. A ^ • 

After normalizing these invariants and repeating the process, two outputs are possible. 
In the first case, the reduced lifted coframe appears to be involutive. Then this coframe is 
the desired set of defining forms for Cont(1Z). In the second case, when the reduced lifted 
coframe does not satisfy Cartan's test, we should use the procedure of prolongation, O 
ch 12]. 

2. Structure and invariants of symmetry groups for nonlinear wave 
equations 

We apply the method described in the previous section to the class of nonlinear wave 
equations (J2J). Denote x 1 = t, x 2 = x, u l = u, u 2 = v, p\ = u t , pi = u x , pf = v t , and 
p\ = v x , The coordinates on 71 are {t, x, u, v, u x , v x }, and the embedding t : 1Z — » J 1 (8) 
is defined by (0). For simplicity in the following computations, we put F(x,u) = 
(a(x,u)b(x,u)Y^ 2 and G(x,u) = (b(x,u)/a(x,u)Y^ 2 , so a(x,u) = F(x,u)/G(x,u) and 
b(x, u) = F(x, u) G(x, u). 

There are three cases to be treated separately: Case A, when F u ^ and G x ^ 0, 
Case B, when G x = 0, and Case C, when F u = 0. 

In the case B system (J2J) has the form u t = F(x, u) (G{u))~ l v x , v t = 
F(x,u)G(u)u x , so the change of variables u = H{u) provided H'{u) = G(u) 
transforms this system into the system u t = F(x, it) v x , v t = F(x, u) u x with F(x, u) = 
F{x,H~ 1 {u)) = F(x,u). Therefore we drop the tildes and conclude that in the case B 
system (J2J) is equivalent to the system 

u t = F(x, u) v x , v t = F(x, u) u x . (7) 
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In the case C system (J2J) has the form u t = F(x) (G(x, w)) -1 v x , v t = 
F(x) G(x,u) u x , so the change of variables x = H(x) provided H'(x) = 1/F(x) 
transforms this system into the system u t = (G(x, u))" 1 v x , v t = G(x,u)u x , with 
G(x,u) = G(H~ l (x),u) = G(x,u). Next, the contact transformation t = v, x = u, 
u = x, and v = t maps the last system to the system in the form uj = F(x, u) v x , 
vj = Fix, u) u x . Thus in the case C system (j2J) is equivalent under a contact 
transformation to the system in the form (J7|) too. 

Let us analyze the system 

u t = F(x, u) (G(x, w)) -1 v x , v t = F(x,u) G(x,u) u x (8) 
in the case A. For brevity we denote 

P = ^W- (9) 

Computing the linear dependence conditions for the reduced forms 6 a , £ l , and erf by 
means of MAPLE, we express the group parameters a\, a\, b\, b\, f xx , f\ 2 , f 2X , /| 2 , g\ x , 
9m 9n, and 9I2 • Particularly, since 

x F{a\-Ga\){a\ + Ga\) det(a$) 2 . 2 - 2 x 

= G(6i-G^)(6i + G^)det(6S) " 2 (m0 ^ >° ' C ' C 
we take a 2 = G~ l a\. Then 

and we take b\ = Fb\. Similarly, we set the coefficients of <j\ at 9 1 , 9 2 , ?, and £ 2 equal 
to zero, and express f xx , f X2 , g\ x , and g\ 2 , respectively. 
Then we obtain 



2 

0~i 



F(a 1 + Gaf)bl 1 , A a i ra t i e 2 2 

(^ + F 6 i )fl i ° 2 (mod^,^,e,e,a 2 



so we take a\ = —Ga\. Now we get 

2Fb? 



a, = 



a\ (mode 1 ,?,?,? 



1_ (b\ + Fb\) 

Since 6 2 7^ (otherwise 6 2 = and det(6*) = 0), we set the coefficient at o\ equal to 1, 
and obtain b\ = — (F^b] + 2 6 2 ). After that, we set the coefficients of a\ at 9 l , 9 2 , ?, 
and £ 2 equal to zero and find f 2X , f 22 , g\ 2 , and g 22 , respectively. This yields 

a\ = 0, aj = o\. (10) 

At the next step we analyze the forms 0^ = t* and ipj = l* reduced by 
setting (Hill and substituting the values of a\, a\, b\, b\, f\ x , f\ 2 , f% x , /| 2 , flrj 2 , 
g 2 2 , and g 22 obtained at the previous step. The form <p\ is semi-basic now, and 
2 = c\ a\ (mod 6* 1 , 6* 2 , ? , £ 2 , a 2 ). So we take c 2 = 0. For the semi-basic form 2 we 
have <p\ = (c 2 + c\) cr^ (mod 1 , 9 2 , ?, £ 2 , erf), so we put c 2 = — c\. Then 

<i _ Fu{P — F Gu x — F v x ) a\ x i^^ r \ oi 02 ( 2\ 
<p2 = AFG 2 a 2 {b\ + Fb 2 ) ^ (^od^0,e), 
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and we take a\ = F u a\(P — F Gu x — Fv x ) F~ l G~ 2 (b{ + Fb 2 2 )- 1 . Then we have the 
semi-basic linear combination ip\ — <p\ + 2 with 

so we take c\ = F Gb\ (P — F Gu x — F v x )~ l (a\)~ l . Similarly, we set the coefficients of 
4>\ and 2 at £ 2 equal to zero, and find f 22 and respectively. Then 

2 _ F 2 (P - FGu x - Fv x ) (P - FGu x + Fv x ) 
01 = AF*G 2 b 2 (b\ + Fb 2 ) * (mod 0,0), 

so we set the coefficient at £ x equal to 1 and find 

x _ Fl (P -FGu x -Fv x )(P-FGu x + F v x ) — AF 4 G 2 (b 2 ) 2 
1 AF 3 G 2 b 2 2 

Then the semi-basic linear combination ip\ + 2 (0 2 — 4>\) gives 
i// + 2^ 2 _ 16F^W^ - - FGu x ) 2 - F 2 v 2 x ) 

^ + 2(^-00 = leP^a} (6|) 2 ^ ( mod ^^^O, 

therefore we put 

, _F tt 3 ((P-FG«,) 2 -F^ 

Go — 



16F 5 G 2 a}(6i) 2 

At the third step, we analyze the reduced structure equations. After absorption, 
we have an essential torsion coefficient at C, 1 A o\ in do\. This coefficient depends on 
/ 2 2 ; we set the coefficient equal to zero and express f 2 2 , while the expression is too long 
to be written in full. Similarly, we express f 2l from the essential torsion coefficient at 
£ 2 A a 2 in da 2 . Then after absorption of torsion in all the structure equations we have 

da\ = Ci A (2 6 l + a 2 ) + ( 2 A (6 1 + a 2 ) + C 3 A 9 2 + Q A (£ 2 + £ 2 ) - £ 2 A a 2 

2F 5 G 2 (b 2 ) 2 (GP x + (Gu x -v x )P u ) q2a 2 



F 2 a\ (P - FGu x + PuJ (P - FGu x - Fv x f 2 ' 

where Ci, C2, C35 an d C4 are 1-forms onTlxH, and the last torsion coefficient is essential. 
There are two possibilities now: P 7^ const and P = const. Denote by V\ the subclass 
of systems (jHJ) such that G x 7^ 0, F u 7^ 0, and P 7^ const. For a system from Pi we set 
the coefficient at 9 2 A a\ in dcr| equal to 1 and obtain 

i = 2P 5 G 2 (6 2 ) 2 (GP a + (G^-^)P M ) 
fll P 2 (P - PGu x + P^) (P - FGu x - Fv x f ' 

Similarly, we set the essential torsion coefficient at 8 1 A 9 2 in d^ 1 equal to zero and find 

2 _P^i 3 : -rai^F^ 
62 " 2P 2 G 

Next, we express g\ 2 from the essential torsion coefficient at 9 1 A £ x in da\. Now the 
essential torsion coefficient at 8 1 A o\ in c?^ 1 has the form 

_ F u (FGP X + PP U ) (P - PGn, + Fv x ) 2 
ix = - 



2P3 (GP x + (Gu x -v x )P u 



>2 
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This function is an invariant of the symmetry group for a system from V\, together with 
its invariant derivatives V^R), i E {1, 6}, defined by dR = V X {R) 9 1 + V 2 (R) 9 2 + 
V 3 (R) (* + V 4 (R) £ 2 + V 5 (R) o~2 + V 6 (R) a\. The invariant V 3 (R) depends on g 2 ^, we 
set T> 3 (R) = and express gf v 

Now all the parameters of the group TC are expressed as functions of x, u, u x , and 
v x . The structure equations of the symmetry group for a system from V\ have the form 

dO l = | (6K 3 + 1 -AK 2 K 3 K 5 -2K 4 ) 9 1 A 6 2 

+\{K A + 1 + 2 K 2 K 3 K 5 - 6 K 3 ) K 2 ~ 1 K 3 1 9 1 A ^ 

+i (AT 4 + 1 + 3 #i# 2 # 3 + 3 # 2 # 3 # 6 - 6 K 3 - 4 K 2 K 3 K 5 ) K 2 ~ l K 3 l 9 l A £ 2 

+ir 2 ir 3 e i A f 7 2 2 -|0 2 Ae 1 +e 2 A^ ) 

de 2 = k 4 e 1 a e 2 - e 1 a f 1 + # 5 # 2 a + # 6 # 2 a «e 2 + # 2 # 3 fl 2 A a\ + (f 1 + £ 2 ) A a 2 , 

df 1 = K X K 2 9 1 A f 1 + \ K X K 2 9 2 A^ + K x ^ A £ 2 , 

d£ 2 = # 2 # 3 9 1 A9 2 -9 1 A^ 1 + [K X K 2 - 1) tf 1 A £ 2 + ± (2 # 3 - 1 + K X K 2 ) 9 2 A £ 2 

+i (# 4 + 1 - 4 # 2 # 3 # 5 - 6 # 3 ) ir 2 - 1 ^ 3 - 1 A e 2 , 
dfji = _if 10 Q 1 /\Q 2 — K u 9 1 A £ 2 + (K A - 1) 9 1 A o\ 

+± (6 K 3 - 1 + 4 # 2 # 3 # 5 - # 4 ) Ko- 1 ^ 1 fl 2 A f 1 

-# 9 fl 2 A £ 2 + |(2 K A - 3 K X K 2 + 4 K 2 K 3 K 5 - 12 K 3 + 2) 9 2 A a\ 

+K 8 e A e - l (AT 4 - 6 K 3 + 1 - # 2 # 3 # 5 ) K~ X K^ ^ A a\ + \ £ A a\ 

+K 7 £ 2 A o\ - K 12 £ 2 A a\ + # 2 # 3 <r£ A 
rfa^ = Km 9 1 A 9 2 + K 17 A 9 1 A + K u 9 1 A £ 2 + (K 4 - K X K 2 + 1) 9 1 A o\ 

+K 18 9 2 A f 1 + #i5 # 2 A £ 2 + § (2 #2# 3 #5 + K 4 - 2) 9 2 A a 2 2 - K m ^ A £ 2 

+K l3 e A a\ - e A a\ - K 20 £ 2 A a 2 , 
where 

K _ 2Fv x (2 F 2 GP + F U G U FP - F UU FGP - G 2 F 2 F XU - F U FGP U + G 2 F X F U F) 

1 (P - u x GF + P Wx ) F U 2 G (-P + + Fv x ) 

= F« (-P + «xGf + F^) (-P + u*GF - -FVc) 2 

2 2P 3 ^ (u x P u G + P X G - P u v x ) 

= v x (PP U + FGP X ) 

3 (FGu x + Fv x - P) (GP X + GP u u x - P u v x ) ' 

while the expressions for K A , K 2Q are too long to be written in full. 

The functions K\, K 2Q are differential invariants of the symmetry group 
Cont(jV) for system (jSJ) from V\. All the other differential invariants of Cont(TZ) are 
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functions of Kj and their invariant derivatives Kjj = Vj(Kj), where for a multi-index 
/ = («i,i 2 , of length #/ = I we denote £>/ = oD i2 o ... oP ip z fe g {1,...,6} 
for k G {1,...,/}. For s > the s i?t order classifying manifold associated with the 
coframe = {9 1 , 9 2 , £ 2 , o^, cr|} and an open subset V in space M 4 with the coordinates 
(x, u, u x , v x ) is 

C {s \6,V) = {(K jJ (x,u,u x ,v x )) | j G {1, ...,20}, #I<s, (i,tt I « l) i),)GV}. (11) 

Since all the functions if^j depend on four variables x, u, u x , and v x , it follows that 
p s = dimC^(0, V) < 4 for all s > 0. Let r = min{s | p s = p s+x = p s+2 = •••} be 
the order of the coframe 6. We have < po < p\ < p 2 < ... < 4. In any case, 
r + 1 < 4. Hence from Theorem 8.19 of ^3] we see that two systems (JEj) from the 
subclass V\ are locally equivalent under a contact transformation if and only if their 
fourth order classifying manifolds fTTf locally overlap. The dimension of Cont(TZ) is 
equal to 6 — dimC^ (0,V). Therefore dimCont(lZ) > 2, as it should be, since every 
system (jHJ) is invariant under the symmetries with infinitesimal generators d/dt and 
d/dv . 

Now we consider the case G x ^ 0, F u ^ 0, and P = m = const. From (jHJ) it follows 
that the system H x = — mF _1 , H u = G is compatible, therefore there exists a function 
H(x, u) such that dH = —mF~ l dx + Gdu. Then the change of variables u = H(x, u), 
v = v — mt maps system (jHj) to the system u t = F(x,u)v x , v t = F(x,u)u x with 
F(x,u) = F(x,u). Dropping tildes, we obtain system 0. Thus in the case P = const 
system (jHJ) is equivalent to system (J7|). 

Let us consider system (J7J). The computations are similar, so we omit them and 
present the results. The structure of the symmetry group for system is different in 
the cases of (\nF) xu ^ and (lnF)^ = 0. We denote by V 2 the subclass of systems 
((Zj) such that (lnF) xu ^ 0. For a system from V 2 all the parameters of the group 7i are 
functions of x, u, u x , and v x . The structure equations for the coframe 6 have the form 

d9 l = (L 3 9 1 +Z 1 )A6 2 + l (3 L 2 L 3 + L 2 - L 4 + L x - L t L 3 ) 9 l A^ + (L 4 9 1 - a\) A £ 2 , 

d9 2 = i L 3 9 1 A 9 2 + 9 1 A £ 2 + § (3 L 2 L 3 - 2 L 2 - L 4 + L x - L X L 3 ) 9 2 A f 1 
-(2 L 2 + Li - L 4 ) # 2 A £ 2 + (f 1 + e 2 ) A cr 2 2 , 

d^ = -l9 1 A^ 1 + 9 2 A^ 1 + L x f 1 A £ 2 , 

d£ 2 = " I 1 A £ 2 + 9 2 A £ 2 + L 2 £} A £ 2 , 

= § (4L 2 + 3L 2 L 3 - L3L1 — 4 L 4 + 6 + 4 Li) (0 1 A 9 2 — £ 2 A a\) + L 7 tf 1 A ^ 
+ (2 £ 2 ^3 ~ y -^2 ~ ^ -^ 2 -^3 + I -^2-^4 + 73 -^2 — I L X L 2 L 3 — 7 £2-^1 — y L 4 + | L X 2 L 3 
+| LxL 4 - § L3L1 + 1 + f Li - \ L 6 - I L\) 9 1 A £ 2 + \ L 3 9 1 A a\ - 2 L 2 9 2 A ^ 
+9 2 A (L 6 e 2 - (L 3 + 1) cr^) + 1 (14 L 2 + 6 L 2 L 3 - 2 Ls^! - 8 L 4 + 3 + 14 L{) £ 2 A o\ 
+L b i l A£ 2 + \ (L 4 - 3L 2 L 3 - 4L 2 - L x + LxL 3 ) ^ A ^ + ^ A a 2 2 , 
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do\ = -| (3 L 2 L 3 + L 2 + | + L x - LxLa - L 4 ) tf 1 A # 2 + i (2 L 2 - 2 L 2 L 3 + 4 L^g 
+2 L 1 L 2 L 3 - 2 LxL 4 - L x + 6 L X L 2 + 6 L 8 + 2 L 2 L 4 + L 4 - 6 - 12 L 2 L 3 + 4 L 2 
+ 12 L 2 L 3 + 2 L 2 ) fl 1 A (e 1 - £ 2 ) + H L 3 + 1) O 1 A a 2 - (L x - L 8 + L 7 - 3 L 2 ) 9 2 A ^ 
+L 8 O 2 A - £3 # 2 A a 2 - I (3 L 2 L 3 + 4 L 2 - L 4 + L x - L^) f 1 A o\ 
+\ (2L 4 - 4L 2 + 1 - A o\ + \ (L 4 - 3L 2 L 3 - 4L 2 - L x + L X L 3 ) £ 2 A ^ 

+ |^5 A e 2 + I (6 L 2 L 3 + 2 2 + 3 + 2L!-2 L X L 3 - 2 L 4 ) £ 2 A a 2 , 
where 

Z/i = (3 v x F u 3 F uu v x F 5 v x F xu F + 5 v x F x F u 3 u x F x F u + 3 F uu u x F 3 u x F u 

+3u x F xu F) {ul-vly l F~ 2 , 
L 2 = (F 2 L lUx + 8 u,F 2 - 8 ^FF m + F 2 !^) F~ 2 (3 « x - 5 <T\ 
^3 = m F u 3 (m 2 - v 2 x ) (6 u 2 x F u L 2 L x - u 2 x F u L\ - 9 m 2 F u L 2 + 8 F x L x v x - 24 F^L 2 , X 

+8 FL x>x v x - 6 F u v 2 x L 2 L x + F u L\v 2 x - 24 F x L 2 v x + 9 i^; 2 ! 2 ) v x 2 {FF xu - F X F U )~ 2 , 
L A = -±F U (u 2 x - v 2 x ) (Qu 2 x F u L 2 2 + 9u 2 x F u L 2 -3u 2 x F u L x - Au 2 x F u L\ + 10 u 2 x F u L 2 L x 

+ 18 ^F X L 2 - 6 u x F x L x - 6 F u v 2 L 2 2 + 10 F x L x v x + 4 F u L\v 2 x + 16 FL^u, 

-9 F u v 2 x L 2 + 3 F u Liv 2 x - 30 F^u, - 10 F u v 2 x L 2 L x ) {FF XU - F X F U )~\ 

while L5, L§ are too long to be written in full. All the differential invariants 
of Cont(TZ) are functions of Lj and their invariant derivatives Ljj = Vj(Lj) = 
T>i x o T>i 2 o ... o X>j ; (Lj) (the operators V>i are not the same as in the case V x \) The 
s th order classifying manifold associated with the coframe and an open subset V is 

C {S \6,V) = {(L jtI (x,u,u x ,v x )) I j G {1, 8}, #/<s, {x,u,u x ,v x )eV}. (12) 

Since all the functions Ljj depend on four variables x, u, u x , and v x , it follows that 
p s = dimC^(6, V) < 4 for all s > 0, and the order r of the coframe 6 satisfies r + 1 < 4 
again. Two systems (J7J) from the subclass V 2 are locally equivalent under a contact 
transformation if and only if their fourth order classifying manifolds (|T2*|) locally overlap, 
and dimCont(TZ) = 6 - dimC (4) (0, V) > 2. 

If (mF) xu = 0, then F(x,u) = S(x) F(u) for arbitrary functions S and F. 
Then the change of variables x = H(x) provided H'(x) = (S^x)) -1 maps the system 
u t = S(x) F{u) v x , v t = S(x) F{u) u x , to the system u t = F(u) v x , v t = F(u) u x . We 
drop the tildes for simplicity of notation and consider the system 

u t = F(u)v x , v t = F(u)u x . (13) 

The computations show that there are three non-equivalent types of systems f|13|) : 
denote by V 3 the subclass of systems flTB^) such that F u 7^ and 
_ _ 4FF 2 F,„ + 4F 2 F. 2 . - 4F 2 F.F,^. - 3F 4 . 
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by P 4 denote the subclass of systems (|T5|) such that F u ^ and Mi = const, finally, by 
V5 denote the subclass of systems (fT3*j) such that F u = 0. 

The subclass V3 is not empty; for example, system (fl~3j) with F{u) = (1 + u 2 )~ l 
belongs to V3. For a system from V3 the structure equations of the symmetry pseudo- 
group after a prolongation have the form 

de l =i) 1 A« 1 -fl 2 A( 1 +( 2 Aff 2 1 1 

do 2 = rj! a e 2 + e 1 a £ 2 - m 2 e 2 a f 1 - (2 m 2 + m 3 ) # 2 a £ 2 + (f 1 + £ 2 ) a of, 

d^M^Af, 
d£ 2 = M 2 A £ 2 , 

rf^ 1 = r/i A o\ + i|2 A ( 2 + Mi 1 A( 1 -2M 2 e 2 A( 1 - M 2 f 1 A ^ + f 1 A a 2 , 
rfa 2 2 = m A a 2 + %A (f 1 + £ 2 ) + Mi 2 A - 2 (M 2 + M 3 ) £ 2 A a 2 , 
^^(Mi-^^A^ 2 , 

dm = Aii A £ 2 + Vi A ^2 + 2 M 2 772 A - V3 A f 1 + (D 4 (Mi) + 2M 2 - M X M 3 ) 9 1 A £ x 
- (2£> 4 (M 2 ) + 4M 2 - Mi) 9 2 A f 1 + (P 4 (M 2 ) + 2Mi - M 2 M 3 - 1) f 1 A 
-(4M 2 + M 3 )^ A a 2 , 

dvs = [i2 A (£* + e 2 ) + »7i A m ~ 3 (M 2 + M 3 ) 773 A £ 2 - (2MiM 2 + 1) 6 2 A £ 2 

+ (4Mi - 2£> 3 (M 2 )(£> 4 (M 1 ) - 1) + 2M 2 (M 2 + M 3 ) - 3) £ 2 A a 2 , 

where 771, ?7 2 , 773, /ii, and /i 2 are 1-forms on TZ x H. The only non-zero reduced 
character, |T3*1 def 11.4], is = 2, therefore the symmetry pseudo-group for system 
(|13|) from "P 3 depends on two arbitrary functions of one variable. The invariants M 2 
and M 3 are defined by M 2 = {2FF uu M hu - FF u M huu - 2F 2 M hu )F- l F- 1 M^ 2 and 
M 3 = — (M 2 D 4 (Mi) + P( 3 ,4) (Mi)), where for an arbitrary function R{u) we have 
ciR = T> 3 (R) + £> 4 (i?)£ 2 with the invariant derivatives D 3 = M{^d/du and £> 4 = 
(1 - 4F 2 M 2 U F- 2 ) M^ We have P 3 (M X ) = 1 and I> 4 (Mi) =' 1 - 4F 2 M 2 U F- 2 . 

Since Mi ^ const, then M 2 and P 4 (Mi) depend on Mi functionally: M 2 = Hi (Mi) 
and P 4 (Mi) = H 2 (Mi). All the other differential invariants can be expressed as 
functions of M 1 . For example, we have £> 3 (M 2 ) = H[{M 1 )V z {M l ) = H[{Mi) and 
P 4 (M 2 ) = P 4 (Mi) V 3 {M 2 ) = F 2 (Mi) Fl(Mi). 

The first order classifying manifold associated with the coframe 6 = 
{9 1 , 9 2 , £ x , £ 2 , <r 2 , of, rji, r}2, rjs} and an open subset C R can be parameterized by 
Mi, M 2 , and £> 4 (Mi) : 

C«(6>,iy) = {(Mi(«),M 2 («),X> 4 (M 1 )(«)) I uG^}. (14) 

Two systems ()13j) from "P 3 are equivalent under a contact transformation iff their 
classifying manifolds ()14j) (locally) overlap, [T3*l Th 15.22], i.e., they have the same 
functions Hi and if 2 . 
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The subclass P 4 is not empty; for example, systems ()13|) with F(u) = 
exp(C arctan(sinh(A u))), F(u) = e", or F{u) = u m , m ^ 0, belong to For a 
system from P4 the structure equations of symmetry pseudo-group after a prolongation 
have the form 

dO 1 =i )l A« 1 -fl 2 A( 1 +( 2 Aff 2 1 , 

de 2 = t) 2 a e 2 - e 1 a e + (e + e 2 ) a a 2 2 , 

^ 1 = (li-1 2 )A(( 1 + 2( 2 ), 
^ 2 = (^-^i)A^ 2 , 

do^ = -2th A (# 2 - ^2) + r]2 A (2 9 2 — a\) + 773 A £ 2 + M 1 9 1 A f 1 + £ x A a 2 , 
da 2 = (2% - ^) A a 2 + 774 A (£ x + £ 2 ) + M x # 2 A £ 2 , 
^^(Mi-^^A^ 2 , 
^^-(Mi-^^A^ 2 , 

dr)3 = /xi A £ 2 + (3 ryi - 2772) A 773 - 2 (M x + 1) - 772) A 9 1 + 4 (771 + 772) A a 2 - 774 A ^ 

+ (3 M x - 4) 6 2 A C 1 + (4 Mi - 3) f 1 A <r\, 
dm = n 2 A (e 1 + e 2 ) - (2^ - 377 2 ) A 774 + (4Mx - 3) £ 2 A a 2 , 

where 771, r/ 2 , 773, 7/4, and /x 2 are 1-forms on 1Z x 7i. The only non-zero reduced 
character is s' x = 2, therefore the symmetry pseudo-group for system (|13J) from V3 
depends on two arbitrary functions of one variable. Since M\ = const, all the other 
differential invariants are equal to zero, and the classifying manifold is a point. Two 
systems from Va are equivalent under a contact transformation iff they have the same 
values of M\. 

A system from P 5 with F(u) = m = const can be transformed to the system 

u t = fx, v t = u x (15) 

by the change of variables t 1— > m~ l t. The symmetry pseudo-group for system ()15|) has 
the structure equations 



d0 l 


= r)i f\9 1 + e Acr^, 


de 2 


= r/2 A 8 2 + e A a 2 2 , 


de 


= %A£ 1 + 77 4 A#\ 


de 


= % A^ + % A# 2 , 


da\ 


= (Vl - A (T2 1 + 777 A 6 1 + 7/ 8 A £\ 


daj 


= (tfc - Vb) Aaj + ri 9 A6 2 + r] 10 A ^ 2 
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where rji, 7710 are 1-forms on 1Z x 7i. The non-zero reduced characters are s[ = 6 and 
s' 2 = 4, therefore the pseudo-group depends on 4 arbitrary functions of two variables. 

The subclasses V3 and V4 are linearizable: the contact transform t = v, x = u, u = 
x, and v — t maps system (|T3j) to the system ui = F(x) v%, v>i = u$. Therefore 

all systems (JBJ) with infinite-dimensional symmetry pseudo-groups are linearizable, cf. 

The results of the computations are summarized in the following 

Theorem : Every system from the class of nonlinear wave equations (0) is equivalent 
under a contact transformation to a system from one of the five invariant subclasses 
Pi, 'Pit P3, P^ an d P§- Pi consists of all systems such that G x 7^ 0, F u 7^ 0, and 
G X F 2 F~ l 7^ const, Vi consists of all systems |7p such that (liaF) xu 7^ 0, P3 consists of 
all systems (Ef such that M x = (4FF 2 F UU + 4F 2 F% U - AF 2 F U F UUU - 3F U 4 ) F~ A ^ const, 
P\ consists of all systems MS\) such that M\ = const, and P§ consists of system Q1B\) . 

Systems from V\ and V2 have finite- dimensional symmetry groups, while systems 
from V3, P4, and P^ are linearizable and have infinite- dimensional symmetry pseudo- 
groups. 

Two systems from one of the subclasses V\, P2, orVj, are equivalent to each other 
under a contact transformation if and only if the classifying manifolds 177]) . JU| ) ; or 
P^p for these systems locally overlap. Two systems from the subclass P 4 are equivalent 
if and only if they have the same constant value of the invariant M\. 

Conclusion 

In this paper, the moving coframe method of [E] is applied to the local equivalence 
problem for a class of systems of nonlinear wave equations under an action of the pseudo- 
group of contact transformations. We have found five invariant subclasses and shown 
that every system of nonlinear wave equations can be transformed to a system from one 
of these subclasses. The structure equations and the differential invariants for all the 
subclasses are found. The solution of the equivalence problem is given in terms of the 
differential invariants. Three of the invariant subclasses consist of linearizable systems 
with infinite-dimensional symmetry pseudo-groups. Therefore all the linearizable cases 
for non-linear wave equations are classified. 
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